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ABSTRACT 

Cappellari (2008) presented a flexible and efficient method to model the stellar kinematics of 
anisotropic axisymmetric and spherical stellar systems. The spherical formalism could be used 
to model the line-of-sight velocity second moments allowing for essentially arbitrary radial 
variations in the anisotropy and general luminous and total density profiles. Here we generalize 
the spherical formalism by providing the expressions for all three components of the projected 
second moments, including the two proper motion components. A reference implementation 
is now included in the public JAM package available at http;//purl.org/cappellari/software. 
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1 INTRODUCTION 

In Cappellari (2008) we used the Jeans (1922) equations to de¬ 
rive the projected second velocity moments for an anisotropic ax¬ 
isymmetric or spherical stellar system for which hoth the luminous 
and total densities are described via the Multi-Gaussian Expan¬ 
sion (MGE, Emsellem, Monnet & Bacon 1994; Cappellari 2002). 
We called the technique the Jeans Anisotropic Modelling (JAM) 
method and provided a reference software implementation^ (in IDL 
and Python). An implementation^ in the C language was provided 
by Watkins et al. (2013). 

In an addendum (Cappellari 2012) we gave explicit expres¬ 
sion for the six projected second moments, including both proper 
motions and radial velocities, for the axisymmetric case (see also 
D’Souza & Rix 2013; Watkins et al. 2013). All projected compo¬ 
nents can be written using a single numerical quadrature and with¬ 
out using special functions. 

In this short note we do the same for the spherical case, namely 
we provide explicit expression for the three components of the pro¬ 
jected second velocity moments. We adopt identical notation and 
coordinates system as in Cappellari (2008), and we refer the reader 
to that paper for details and definitions. 


2 JEANS SOLUTION WITH PROPER MOTIONS 

We assume spherical symmetry and constant anisotropy for each 
individual MGE component. The Jeans equation can then be written 
(Binney & Tremaine 2008, equation 4.215) 
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where Vq = for symmetry and we defined = 1 — Vg/vr. 

We use equation (19) of van der Marel & Anderson (2010) 
which provide the projection expressions for all three components 
of the velocity second moments, including the proper motions (see 
also Strigari, Bullock & Kaplinghat 2007). We then follow the same 
steps and definitions as Cappellari (2008, section 3.2.1) to write all 
three projected second velocity moments as follows 


T.vUR) = 2G 





dr, 


( 2 ) 


where (i) a — los for the line-of-sight velocity (ii) a = pmr for 
the radial proper motion, measured from the projected centre of the 
system, and (iii) a = pmt for the tangential proper motion, respec¬ 
tively and we defined 

QUr) = 1 - P (R/rf (3) 

Qpmrir) = 1 - P + P {R/rf (4) 

Qpmt(r) = 1-/3. (5) 


Integrating by parts one of the two integrals disappears and 
all three projected second moments can still be written as in equa¬ 
tion (42) of Cappellari (2008) 

T.'^iR) = G j (6) 

As shown in Section 3, when using the MGE parametrization, the 
evaluation of this expression requires a single numerical quadrature 
and some special functions. Eor the line-of-sight component the ex¬ 
pression for was given by equation (43) of Cappellari (2008) 
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(see also Mamon & Lokas 2005), where F is the Gamma function 
(Abramowitz & Stegun 1964, equation 6.1.1) and is the incom¬ 
plete Beta function (Abramowitz & Stegun 1964, equation 6.6.1), 
for which efficient routines exist in virtually any language. The cor¬ 
responding expressions to use for the radial and tangential proper 
motion components are 
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Specific expressions can be obtained for j3 = ±1/2, where the 
Bm function is divergent, but in real applications it is sufficient to 
perturb /3 by an insignificant amount to avoid the singularity. In the 
isotropic limit all three components become equal 


lim ,^ios — lim J ~pmr — 11m R pmt — 2 
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and Equation (6) reduces to equation (29) of Tremaine et al. (1994). 



Figure 1. MGE fit to a de Vaucouleurs (1948) profile. We used 20 
Gaussians and the MGE_FIT_lD routine to describe the profile in the range 
— 3 < logiJ/i?B < 2 with a maximum relative error of about 0.2%. The 
resulting MGE model contains 99.998% of the analytic total mass. 


3 MGE SPHERICAL JEANS SOLUTION 


Here we apply the general spherical Jeans solution with constant 
anisotropy to derive an actual solution for a stellar system in which 
both the luminous density and the total one are described by the 
MGE parametrization. In this case the surface brightness E*, the 
luminosity density zzj. and the total density pj for each individual 
Gaussian are given by (Bendinelli 1991) 
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The mass of a Gaussian contained within the spherical radius r is 
given by equation (49) of Cappellari (2008) 


Mj{r) = M, 




(14) 


with erf(a:) the error function (Abramowitz & Stegun 1964, equa¬ 
tion 7.1.1). 

The projected second velocity moments for the whole MGE 
model, summed over all the N luminous and M massive Gaussians, 
for any of the three velocity second moment components, are still 
given by equation (50) of Cappellari (2008) 


T.vI{R)^gI 


M 


f=i 


dr, (15) 


where Vkir) is given by Equation (12), Mj{r) is given by Equa¬ 
tion (14), and J-a,k is obtained by replacing the /3 parameter in 
Equation (7)-(9) with the anisotropy (3k of each luminous Gaussian 
component of the MGE. 

The formalism presented in this section was implemented in 
an updated version of the public JAM package (see footnote 1). 


4 APPLICATION 

As an illustration of the general behaviour of the second moments, 
we used the Python version of the MGE_FIT_lD one-dimensional 
fitting routine in the MGE_FIT_SECTORS package (Cappellari 2002, 
see footnote 1) to obtain an accurate MGE description of a de Vau¬ 
couleurs (1948) R^^'^ profile (Figure 1). We then used the spherical 
JAM_SPH_RMS routine in the JAM package (see footnote 1) to calcu¬ 
late the predicted velocity second moments for different (here con¬ 
stant) anisotropies, using Equation (15). The model assumes that 
mass follows light except for the presence of a central black hole 
with mass M, = 0.005 x M, where M is the total mass of the 
system. The fractional value is the typical observed one from Kor- 
mendy & Ho (2013). The resulting predicted second velocity mo¬ 
ment profiles are shown in the top panel of Figure 2. The bottom 
panel is like the top one, for a ‘cored’ profile u{r) oc (1 ± r)~‘^ 
with asymptotically constant luminosity density at small radii. 

For the line-of-sight component. Figure 2 illustrates the well 
known mass-anisotropy degeneracy: the profile changes signifi¬ 
cantly at fixed M,. If the anisotropy is unknown, M, cannot be 
measured from these spherical models. However, using proper mo¬ 
tions one can constrain the anisotropy directly and consequently 
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Figure 2. Projected velocity second moments v‘^(R) as a function of ra¬ 
dius, normalized by the half-light radius Rc, for a range of anisotropies 
0, for spherical models with a central black hole mass M* of 0.5% of 
the total mass. The top panel is for a de Vaucouleurs (1948) surface 
brightness profile, while the bottom panel is for a ‘cored’ luminosity den¬ 
sity profile, with asymptotically constant nuclear density. The three adopted 
anisotropies (3 were chosen to correspond to axial ratios of the velocity el¬ 
lipsoid of = 3/4, 1 and 4/3 respectively. Here a. = los, pmr 

or pmt for the line-of-sight, the radial or tangential proper motions respec¬ 
tively. The presence of a core enhances the observable effects of /3 variations 
and reduces the ability of the proper motions to measure /3, especially for 
/3 < 0. This makes the determination of mass profiles more challenging in 
cored than in cusped luminosity profiles. 





measure M, using the second velocity moments alone (van der 
Marel & Anderson 2010). The bottom panel also shows that mod¬ 
els are more challenging for more shallow inner profiles due to the 
reduced sensitivity of the ratio to /3 variations. 

In Figure 2 all three components are plotted in the same units 
as for the radial velocity. In a real application of the method, the 
velocities of the proper motion components need to he converted 
into proper motions units, using the distance of the system. The 
distance is a free model parameter and this allows one to measure 
kinematical distances (e.g. van de Ven et al. 2006; van den Bosch 
et al. 2006; van der Marel & Anderson 2010). 

By assigning different anisotropies Pk to the different Gaus- 
sians one can describe essentially arbitrary smooth radial variations 
of the anisotropy (e.g. Cappellari et al. 2009, section 3.3). While, by 
using different Gaussians for the luminous and total density, one can 
describe general dark halo profiles (e.g. Cappellari et al. 2015). In 


practice, for spherical geometry, one can perform a high-accuracy 
MGE fit to some parametric description of the dark and luminous 
profiles, using the MGE_FIT_lD fitting program of Cappellari (2002, 
see footnote 1) as shown in Figure 1. 


5 CONCLUSIONS 

Proper motion data have the fundamental advantage over the line- 
of-sight quantities that they allow one to break, using the velocity 
second moments alone, the mass-anisotropy degeneracy affecting 
the recovery of mass profiles in spherical systems (Binney & Ma- 
mon 1982). The proposed formalism is becoming useful to model 
the growing number of proper motion data becoming available. 
These are now being provided mainly by the Hubble Space Tele¬ 
scope (e.g. Watkins et al. 2015). In the near future proper motion 
data for the Milky Way satellites will also be provided by the GAIA 
spacecraft, while in the more distant future a major step forward 
in proper motion determinations will be made by EUCLID space 
mission. 


ACKNOWLEDGEMENTS 

I acknowledge support from a Royal Society University Research 
Fellowship. This paper made use of matplotlib (Hunter 2007). 


REFERENCES 

Abramowitz M., Stegun I. A., 1964, Handbook of mathematical 
functions (Reprinted 1972), Dover Books on Advanced Mathe¬ 
matics. Dover, New York 
Bendinelli O., 1991, ApJ, 366, 599 
Binney J., Mamon G. A., 1982, MNRAS, 200, 361 
Binney J., Tremaine S., 2008, Galactic Dynamics: Second Edition. 

Princeton University Press, Princeton, NI 
Cappellari M., 2002, MNRAS, 333, 400 
Cappellari M., 2008, MNRAS, 390, 71 
Cappellari M., 2012, e-print (arXiv:1211.7009) 

Cappellari M., Neumayer N., Reunanen I., van der Werf P. P, de 
Zeeuw P. T, Rix H.-W., 2009, MNRAS, 394, 660 
Cappellari M. et al., 2015, ApJL, in press (arXiv: 1504.00075) 
de Vaucouleurs G., 1948, Annales d’Astrophysique, 11, 247 
D’Souza R., Rix H.-W., 2013, MNRAS, 429, 1887 
Emsellem E., Monnet G., Bacon R., 1994, A&A, 285, 723 
Hunter J. D., 2007, Computing In Science & Engineering, 9, 90 
Jeans J. H., 1922, MNRAS, 82, 122 
Kormendy J., Ho L. C., 2013, ARA&A, 51, 511 
Mamon G. A., Lokas E. L., 2005, MNRAS, 363, 705 
Strigari L. E., Bullock J. S., Kaplinghat M., 2007, ApJ, 657, LI 
Tremaine S., Richstone D. O., Byun Y.-I., Dressier A., Faber 
S. M., Grillmair C., Kormendy J., Lauer T. R., 1994, AJ, 107, 
634 

van de Ven G., van den Bosch R. C. E., Verolme E. K., de Zeeuw 
P. T, 2006, A&A, 445, 513 

van den Bosch R., de Zeeuw T, Gebhardt K., Noyola E., van de 
Ven G., 2006, ApJ, 641, 852 
van der Marel R. P, Anderson J., 2010, ApJ, 710, 1063 
Watkins L. L., van de Ven G., den Brok M., van den Bosch 
R. C. E., 2013, MNRAS, 436, 2598 
Watkins L. L., van der Marel R. P, Bellini A., Anderson J., 2015, 
ApJ, 803, 29 


© 2015 RAS, MNRAS 000, 1-3 


















